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Abstract. General relativity predicts that energy and momentum conservation 
laws hold and that preferred frames do not exist. The parametrised post-Newtonian 
formalism (PPN) phenomcnologically quantifies possible deviations from general 
relativity. The PPN parameter 03 (which identically vanishes in general relativity) 
plays a dual role in that it is associated both with a violation of the momentum 
conservation law, and with the existence of a preferred frame. By considering the 
effects of 03 7^ in certain binary pulsar systems, it is shown that 0:3 < 2.2 X 10 — 20 
(90% CL). This limit improves on previous results by several orders of magnitude, 
and shows that pulsar tests of 03 rank (together with Hughes-Drever-type tests of 
local Lorentz invariance) among the most precise null experiments of physics. 



1. Introduction 

Conservation laws have long been a foundation of physics and careful tests of their 
validity have led to important discoveries, such as neutrinos. Most theories of gravity, 
starting with general relativity, are based on an invariant action principle and therefore 
predict that the energy and momentum of isolated gravitating systems are conserved. 
However, direct experimental constraints on those conservation laws are far less 
extensive than for the other fundamental forces of nature. Another cornerstone 
of physics is the absence of preferred frames in local experiments (Local Lorentz 
Invariance). This property may be violated in certain theories of gravity allowing 
for the existence of long-range vector fields. 

The parametrised post-Newtonian formalism (PPN) phcnomenologically quantifies 
possible deviations from general relativity. Such deviations are measured by a set of 
ten parameters^ — — 1, £, ax, a.2, a^, £i, (2, C3> C4 0- To eacn °f these parameters 
is associated a class of observable, non-Einsteinian effects. The parameter 0:3 plays a 
special role in that it is associated both with a violation of the momentum conservation 
law, and with the existence of preferred frames («3 = in general relativity) . 

The most striking observable effect induced by a putative 0*3 ^ is the existence 
of a non-zero self-acceleration for a rotating body, perpendicular to its spin axis and 
absolute velocity j|, g, Qj. As a consequence, the total acceleration of a body, which 
is a member of a gravitating system, can be decomposed as 



a — a Newt + a nbody + a a 3 



(1) 
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where aN CWt is a two-body Newtonian-like 1/R 2 force (modified by possible 
equivalence-principle- violation effects), a n body denotes relativistic n-body effects, 
and where a Q3 is the self-acceleration resulting from the violation of momentum 
conservation associated with 013 7^ J^, ||, For a nearly spherical body, rotating 
with angular velocity f2, 

a Aa 3 = --"a^-o (w + v) x n (2) 
d rriAC 

where E A iav is the gravitational self energy of body A having mass tua, moving with 
velocity w + v, with respect to the absolute rest frame. For later convenience, we 
have decomposed this velocity into the absolute velocity w of the centre of mass 
of the considered system, and the peculiar velocity v of body A with respect to 
the centre of mass frame. Instead of the gravitational self energy Ea, one can 
introduce the (dimensionless) "compactness" parameter ca = —2d\nrriA/d\nG which 
is approximatively given by ca = —2E A Tav /mAC 2 ||. We see from equation @j that, 
given an absolute velocity w, FIca = 2wca/P defines a figure of merit for the selection 
of bodies testing 013. For the sun ca ~ 10~ 6 , and for a white dwarf ca ~ 10 -4 . By 
contrast, neutron stars have compactness parameters of order unity. Considering a 
range of equations of state for neutron stars, Ref. |Q found that ca takes a median 
value of O.21m^/M0. Therefore, rapidly rotating neutron stars are, by a very large 
factor, the best objects for constraining 013. In the present paper, we shall show that 
existing data on certain long-period, quasi-circular binary pulsar systems allow one to 
constrain 03 at the 10~ 20 level. 



2. Previous Limits 



As shown by Nordtvedt and Will [g, [j]] , a non-zero a.3 induces a contribution to the 
perihelion precession of the planets in the solar system. The two planets with the 
best measurements of periastron advance were Earth and Mercury. By combining the 
observations for two planets it is possible to eliminate the terms involving the well- 
known Eddington parameters 7 and 0, obtaining |49«i — 012 — 6.3 x 10 5 a3 — 2.2£| < 0.1 
0. Using the limits on other PPN parameters a limit of 1 0:3 1 < 2 x 1CP 7 was thus 
obtained. 

A tighter limit on has been obtained by considering the effect of the acceleration 
(||) on the observed pulse periods of isolated pulsars. The observed pulse period 
P = Pq(1 + v r /c), contains contributions from both the intrinsic pulse period Pq 
and the Dopplcr effect resulting from the radial velocity v r . Consequently any radial 
acceleration a r contributes to the observed period derivative P, 

t = + (3 ) 
P P c K ' 

Self accelerations are directed perpendicular to both the spin axis and the absolute 
velocity of the spinning body. If self accelerations were contributing strongly to the 
observed period derivatives of pulsars, roughly equal numbers of positive and negative 
observed period derivatives would be expected, since the spin axes and therefore the 
self accelerations are randomly oriented. However the observed distribution (excluding 
those pulsars in globular clusters) contains only positive period derivatives, allowing 
a limit to be placed on 013. This was done initially using normal pulsars EJ to obtain 
a limit of | CK3 1 < 2 x 10~ 10 and more recently using millisecond pulsars H to obtain 
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a limit of |as| < 5 x 10 -16 . The later case incorrectly included binary pulsars in the 
sample. [As will be shown in the following section the effects of self accelerations in a 
binary system are more complicated.] Nevertheless, if one restricts oneself to a sample 
of isolated millisecond pulsars, a limit quantitatively similar to that above, though 
slightly less stringent, can be obtained. 



3. Polarisation of Binary Orbits 

When considering a binary system there is a perturbing self acceleration felt by each 
body. These self accelerations perturb both the centre of mass motion of the system, 
and the relative orbital dynamics. We focus on the perturbations of the relative 
motion which turn out to be a much more sensitive probe of a^. The perturbation 
of the relative acceleration a Q3 = a Aoi3 — &Ba 3 , where A labels the pulsar and B its 
companion, can be written as 

a^a^ + a*™ (4) 
where 

a Stark = ^ W X (c A n A - C B n B ) , (5) 



a 



Zeeman 
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v x (x B c A ft A + x A c B fl B ) (6) 



Here c Al £l A , x A — m A /(m A + m B ) are the compactness, angular velocity and mass 
fraction for the pulsar. The names given to the perturbing accelerations (^|) and (^|) 
have been chosen by analogy with the well-known effects of constant external electric 
and magnetic fields on the classical dynamics of an atom. These two perturbations 
have very distinct effects on the relative orbital dynamics. 

By a generalisation of Larmor's theorem, the "Zeeman" contribution @ is easily 
seen to cause a slow overall precession of the orbit with angular velocity 

^precession = ~ — {x B C A Q. A + X A C B fl B ) . (7) 

The theory of the formation of binary pulsars leads us to expect that, to a first 
approximation, the spin vectors £l A , £l B are parallel to the orbital angular momentum. 
This means that the effects of the perturbation (Q) give rise to non-observable, very 
small additional contributions to the orbital period and the periastron precession. 

By contrast, the "Stark" contribution which represents a constant perturbing 
acceleration, leads to a forced eccentricity, polarising the orbit along a fixed direction. 
Such a "Stark" polarisation (which is the DC analogue of the Nordtvedt polarisation of 
the lunar orbit pj ) has already been studied in certain binary pulsars when considering 
the effect of a violation of the strong equivalence principle 0. The best systems for 
constraining the polarisation effects caused by (|^) are the long-orbital-period, quasi- 
circular binary millisecond pulsars with white dwarf companions. For a neutron star 
and white dwarf c B <ti c A and £l B <C 0, A , so that equation (|J) reduces to 

a Sf k ~ ^- CA w x n A . (8) 

For this acceleration to be in a fixed direction, w and tt A must remain fixed 
in space. Any precession of the pulsar spin axes will be very slow, if it occurs 
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at all. The absolute velocity of pulsars in our galactic neighbourhood is taken to 
be |w| = 369 km s _1 , as determined from observation of the cosmic microwave 
background j|] . Some recent results Q question the validity of the cosmic microwave 
background as an absolute reference frame. As a sample of pulsars are used, only the 
magnitude of w is important and this is similar in both cases. 

It is also required that the centre of mass of the binary system does not move 
appreciably in the Galaxy during the build up of the polarisation induced by (||), 
otherwise w will vary and the polarisation force (^|) must be decomposed in a sum of 
monochromatic terms. This is equivalent to a corresponding requirement for the test 
of the strong equivalence principle Q] , ie 

^ = ( Pb ) (^-) « 1 (9) 
ojr Vl364(days)y \l.7M Q ) y ' 

where WQai is the angular rotation rate of the Galaxy near the sun, w R is the relativistic 
periastron advance and M — tua + rriB- 

The equation of motion for a binary system with the perturbation (5) is very 
similar to the equation of motion for the strong equivalence principle test [7 

d 2 r GMv stark 

^2 +^= R + a "s ark (10) 

where R contains the relativistic forces responsible for periastron precession. We 
can use the results of [Q. Considering systems for which e <C 1 while letting 
w x CI a = wQa sin/3e„ (where e„ is a unit vector) and Qa = 2t:/P gives the solution 



e(t) =e R (t)+e a3 = e R (t) + "s^-^-^g sm/3e u . (11) 

Here e(t) is the total, observable eccentricity vector, e R (t) is the intrinsic eccentricity 
vector, which rotates with angular frequency tu R due to relativistic periastron 
precession and e Q3 is the forced eccentricity vector caused by a% ^ 0. It is assumed 
that since considerable mass transfer has taken place, the orbital and spin angular 
momenta are aligned, which results in the forcing term (^|) due to 03 being parallel 
to the orbital plane. Given an observed eccentricity e b s = l e (0l measured at some 
time t, equation (^) may be used to place a limit on 0:3. From equation ( |Tl| ) it can 
be seen that 

PljeP (12) 

defines a figure of merit for selecting which binary pulsar systems will provide the best 
limit. This figure of merit differs from the one (P? /e) selecting the best tests of the 
strong equivalence principle. 



4. Confidence level of the limit 

Unfortunately, in equation (|ll]) there are two parameters which are not constrained 
by the observations, (3 and e R (t). Firstly, the angle f3 may be very small causing 03 
effects to contribute weakly to the observed eccentricity. Secondly, there is a finite 
probability that e R (t) might have cancelled most of the term due to leaving only a 
small observed eccentricity. In order to impose a limit on 013 at a particular confidence 
level a quantitative assessment of the probability of these independent effects is needed. 
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The problem of en(t) potentially cancelling a second term was first considered 
by Damour and Schafer who decomposed the problem into three cases. When 
e n(t) 3> e a 3 or en(t) <§; e Q3 it is sufficient to use the observed eccentricity e bs for 
obtaining a limit. When e^(t) ~ e as and the angle between e^(i) and the opposite of 
e Q3 is 0, then e bs = ea 3 2| sin(0/2)| was used. A slightly more precise approach was 
taken by Wex |l(| who defined a function 

|sin0| if \6\ < tt/2 
1 if |0| > tt/2 



5(0) 



(13) 



Here the angle between en(t) and — e a;! is taken to vary between — 7r and n. Given 



the magnitude of the forced eccentricity e c 



and the value of the angle 0, it 



is easily seen that when the magnitude of the intrinsic eccentricity \en\ is allowed to 
vary one can write the inequality 

e obs > e a3 5(0) . (14) 

[Indeed, when |0| < tt/2 the worst cancellation arises when the resulting eccentricity 
vector e = e^f + e Q3 is perpendicular to e^, while when |0| > tt/2 the worst case is 
|ejj j = 0.] Using the inequality Jli| ) and rearranging equation ( pl| ) gives the upper 
bound 

24^ eP GM 1 = Kj 
la3l -c A wP£ c 2 5(0) sin 13 ~ 5(0) sin (3 1 ' 

Given the quantity Ki determined from observed parameters, the probability that \a^\ 
is greater than some number a is smaller than 

F{Ki/a) = Prob[5(0) sin /3 < Kja] . (16) 

For sufficiently old systems, eu(t) has had the time to make many turns and we can 
consider that the orbital phase is randomly distributed over the interval [—w, n] ■ [The 
pulsars included in our sample below have old white dwarf companions, confirming 
that this oldness condition is satisfied.] Given w, the spin vector Q,a can point in an 
arbitrary direction on the unit sphere. Therefore the polar angle (3 is distributed with 
the probability law ^ sin (3dj3. In other words, the quantity cp — cos (3 is distributed 
uniformly over [—1,1]. Hence, since the distributions are symmetric about = and 
c p = 0, 

|0|<7T r T ^ -\ O<0<7T 



F 



K, 



Prob 



c} < 



4Prob 



|Cjs|<l 



c 2 < 



(17) 



0<Cf)<l 



F 



Ki 



= 1-1/1- 



Vi-(^/«) 2 



dcp arcsin 



Kj/a 



(18) 



Probabilities evaluated from this function using the Mathematica package to 
numerically integrate the second term are shown in Table |l| for a range of values 
of Ki/a. 

Since the observational results for pulsars are independent, the combined 
probability, given the existence of several pulsars with small values for Ki, satisfies 
the inequality 



1 - (CL/100) = Prob[|a 3 | > a] <~[[F{K t /c 



(19) 



() 



Table 1. The probability that S(0) sin/3 < Ki/a 



IS 1 — 

J\i/ a 


Probability 


0.00 


0.000 


0.10 


0.045 


0.20 


0.092 


0.30 


0.143 


0.40 


0.201 


0.50 


0.267 


0.60 


0.343 


0.70 


0.432 


0.80 


0.540 


0.85 


0.604 


0.90 


0.681 


0.93 


0.736 


0.96 


0.804 


0.97 


0.831 


0.98 


0.864 


0.99 


0.905 


1.00 


1.000 



Here, we have introduced the confidence level (CL) with which one can reject the 
hypothesis that \a^\ is larger than some given a. 

When the product on the right-hand side is equal to 0.1, |ct3 1 < a at better than 
the 90% confidence level. Note that when Ki > a, F{Ki/a) = 1, so that including 
such a pulsar does not improve the limit. 



5. The Limit 

In the most recent compilation of fast pulsars in binary systems there are 

many systems which satisfy the constraints discussed in the previous two sections. 
Using the result of || for the median value of ca quoted above, one finds Ki = 
1.925 x 10~ 13 (l+rriB /niA)eP/ P£ where the pulsar period P is expressed in milliseconds 
and the orbital period P& in days. After evaluating Ki for each of those binary 
pulsars, the objects with the smallest Ki were chosen. They are listed in Table g. 
The probability that S(9) sin (3 < Ki/a was evaluated using equation ( p^|) for several 
values of a which are shown in Table ||. The combined probabilities and confidence 
levels from using several pulsars were then determined using equation © (See Table 

D- 

Using, for instance, the 90% confidence level gives a limit \a^\ < 2.2 x 10 20 . This 
is more than a factor of 10 4 better than the previous limit Q , and by far the tightest 
limit on any of the PPN parameters. 

The only other ultra-high-precision null experiments giving limits of order 10 -20 
on a dimensionless theoretical parameter we are aware of are some recent Hughes- 
Drever-type tests of local Lorentz invariance |l2|, [l3|, Q. Figure 14.2 of |l[] shows 
the limits on the parameter S = Cq/c 2 — 1, where cq is the limiting speed of massive 
particles, and c e the speed of light. It is remarkable that tests involving binary pulsars 
can rank, together with modern laser-cooled trapped atom experiments, among the 
most precise null experiments of physics. 
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Table 2. F(Ki/a) for the test systems. Ki and a are in units of 10 



Pulsar 


Ki 


a = 3.7 


a = 2.9 


a = 2.4 


a = 2.2 


a = 1.95 


B1855+09 


17.6 


1.00 


1.00 


1.00 


1.00 


1.00 


J2317+1439 


12.6 


1.00 


1.00 


1.00 


1.00 


1.00 


J1455-3330 


5.45 


1.00 


1.00 


1.00 


1.00 


1.00 


B1953+29 


3.24 


0.64 


1.00 


1.00 


1.00 


1.00 


J1643-1224 


2.30 


0.36 


0.54 


0.80 


1.00 


1.00 


J1640+2224 


1.91 


0.28 


0.40 


0.54 


0.62 


0.86 


J2229+2643 


1.83 


0.27 


0.37 


0.50 


0.58 


0.76 


J2019+2425 


1.79 


0.26 


0.36 


0.48 


0.55 


0.71 


J1713+0747 


1.76 


0.25 


0.35 


0.47 


0.54 


0.68 






0.001 


0.01 


0.05 


0.10 


0.32 


CL > 




99.9% 


99% 


95% 


90% 


68% 
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